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ABSTRACT 

The superb spatial resolution of Chandra has allowed us to detect a 20"-long tail behind the Geminga 
pulsar, with a hard spectrum (photon index F = 1.0 ±0.2) and a luminosity (1.3 ± 0.2) x 10 29 ergs s _1 in 
the 0.5-8 keV band, for an assumed distance of 200 pc. The tail could be either a pulsar jet, confined by 
a toroidal magnetic field of ~ 100 /iG, or it can be associated with the shocked relativistic wind behind 
the supersonically moving pulsar confined by the ram pressure of the oncoming interstellar medium. We 
also detected an arc-like structure 5"-7" ahead of the pulsar, extended perpendicular to the tail, with 
a factor of 3 lower luminosity. We see a 3 a enhancement in the Chandra image apparently connecting 
the arc with the southern outer tail that has been possibly detected with XMM-Newton. The observed 
structures imply that the Geminga's pulsar wind is intrinsically anisotropic. 

Subject headings: pulsars: individual (Geminga) - stars: neutron — stars: winds, outflows 



1. INTRODUCTION 

It has been recognized long ago that pulsars steadily 
lose their rotational energy via relativistic pulsar winds 
(PWs). The PW shocks in the ambient medium and forms 
a pulsar-wind nebula (PWN) which emits synchrotron ra- 
diation (Rees & Gunn 1974). An isotropic outflow from a 
stationary pulsar in a uniform medium forms a spherical 
termination shock (TS) at a radius R s ~ (ij/47rcp am b) 1 / 2 , 
where E is the pulsar's spin-down power and p am b is the 
ambient pressure. For a pulsar moving at a speed v with 
respect to the medium, the TS shape is distorted by the 
ram pressure, p ram = Pamb^ 2 . At supersonic speeds, when 
the ram pressure exceeds the ambient pressure, the TS of 
an isotropic PW acquires a bullet-like shape (Bucciantini, 
Amato, & Del Zanna 2005 [hereafter B05], and references 
therein) with a distance i?h — (i?/47rcp ram ) 1 ^ 2 between the 
pulsar and the bullet head. The shocked PW is confined 
between the TS and the contact discontinuity (CD) surface 
which separates the shocked PW from the shocked ambi- 
ent medium between the CD and the forward bow shock 
(FBS). The shape of the shocks and the overall appearance 
of the PWN depend on the interplay of E, p am b, p am b, and 
v. Thus, PWNe produced by fast-moving pulsars provide 
an important diagnostic tool for studying PWs, the ambi- 
ent medium, and pulsar velocities. 

X-ray PWNe have been observed around about 30 pul- 
sars. High-resolution observations with Chandra show 
that the PWN structure is never a simple spherical shell. 
In particular, young PWNe (e.g., the Crab and Vela 
PWNe; Weisskopf et al. 2000; Pavlov et al. 2003) are 
often approximately axially symmetric, with jets along 
the symmetry axis (which apparently coincides with the 
pulsar's spin axis) and a torus-like structure around the 
axis (an equatorial PW outflow). There are some X-ray 
PWNe whose cometary shape is clearly caused by the pul- 
sar motion. Some of them (e.g., the Mouse PWN, pow- 
ered by the young PSR J1747-2958; Gaensler et al. 2004) 



are confined within a bow-like boundary, but they do not 
show a shell structure, in contrast to the sharp H Q "bows" 
emitted by the shocked ambient gas at FBSs produced by 
some pulsars. Chandra observations of the pulsars B1757- 
24 (Kaspi et al. 2001), B1957±20 (Stappers et al. 2003), 
J1509-5859, and J1809-1917 (Sanwal et al. 2005) have 
revealed elongated structures (of lengths ~ 0.05-0.5 pc) 
which look like "trails" behind moving pulsars. An ex- 
ceptionally long (~ 1.6 pc) trail was found behind PSR 
B1929±10 by Wang, Li, & Begelman (1993). The origin 
of these trails has not been firmly established (see §4). 

An intriguing PWN structure has been recently reported 
by Caraveo et al. (2003; hereafter C03) from an XMM- 
Newton observation of the middle-aged pulsar Geminga 
(t = 340 kyr, E = 3.3 x 10 34 ergs s" 1 ). They found two 
2'-long tails behind the pulsar, with a luminosity of ~ 10 29 
ergs s _1 in the 0.3-5 keV band. They suggested that these 
are the tails of a bow shock generated by pulsar's motion 
and predicted the head of the bow-shock to be at an angu- 
lar distance of 20"-30" ahead of the pulsar. In this paper, 
we present the results of our high-resolution Chandra ob- 
servation of the Geminga PWN 3 and confront them with 
the results of the XMM-Newton observation and theoret- 
ical models. 



2. OBSERVATIONS AND DATA REDUCTION 

Chandra observed Geminga on 2004 February 7 for 
19.9 ks (18.8 ks good exposure time). The pulsar was 
positioned on ACIS-S3 chip with a standard Y offset of 
— 0.'33. A 1/8 subarray mode (data taken from a V x 8' 
region of each CCD, with a frame time of 0.7 seconds) was 
used to mitigate the pilcup in the pulsar image. We used 
the Chandra Interactive Analysis of Observations (CIA0) 
software (ver.3.2.2; CALDB ver.3.1.0) for the analysis, start- 
ing from the level 1 event files to correct for charge trans- 
fer inefficiency. We applied the standard grade filtering 
and removed events with energies > 8 keV to reduce the 
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Fig. 1. — Left: XMM-Newton EPIC-PN image of a 4' X 4' region around Geminga, smoothed with a 15" FWHM Gaussian. The dashed 
curves show two bow-shock models by Wilkin (1996), with stand-off distances 20" and 7"; the former fits the shape of the two outer tails, 
while the latter corresponds to the actually observed arc ahead of the moving pulsar (see the right panel). The arrow shows the direction 
of pulsar's proper motion. The black circle indicates the position of the point source X resolved in the Chandra image (see Fig. 2). Right: 
Chandra ACIS-S3 image of a 1' X 1' region around Geminga, showing an extended emission at 5"— 7" ahead of the pulsar (the arc) and a 
20"-long axial tail. 



1 x 10~ 6 counts s _1 arcsec~ 2 in the 0.3-8 keV band. Spec- 
tral analysis of the tail emission is hindered by the strong, 
nonuniform background. Using the absorbed power-law 
spectral model with the parameters estimated by C03 from 
the MOS data (photon index T = 1.6, hydrogen column 
density iV H = 1.1 x 10 20 cm" 2 ), we obtain the surface 
brightness (intensity) of the outer tails /outer ~3x 10 -18 
2 • "' ' - ! and the total unabsorbed flux 



ergs cm s arcsec 
Fouta, ~ 2 x 10- 14 



ergs cm 2 s 1 . 
1 x 10 29 rf2oo ergss -1 



This gives the lumi- 
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background. Since no events with E < 0.45 keV were 
telemetered from the S3 chip, we use the 0.5-8 keV band 
for further analysis. To account for nonuniform exposure 
and nonuniform CCD response near the node boundaries, 
we applied the exposure map correction, using the CIA0 
script merge_a.ll. 

We also reanalyzed the XMM-Newton EPIC data taken 
on 2002 April 4-5 (103 ks exposure, with EPIC-PN in 
Small Window mode and EPIC-MOS in Prime Full Win- 
dow mode). We used SAS v. 6.1 to reprocess and analyze 
the data. The effective exposure times, after removing the 
periods of high background, are about 68 ks for EPIC-PN 
and 82 ks for MOS1 and MOS2. A detailed analysis of the 
MOS data has been presented by C03. We will supplement 
those results with the analysis of the EPIC-PN data. 

3. OBSERVATIONAL RESULTS 

3.1. XMM-Newton EPIC-PN Results 

Smoothed MOS and PN images show two elongated 
patchy structures, stretched in the direction opposite to 
the pulsar's proper motion. We will call them outer tails 
to distinguish from the short inner tail discovered with 
Chandra (Sanwal et al. 2004; see §3.2). Figure 1 (left 
panel) shows these structures, seen up to about 3' from the 
pulsar (~ 0.17<i2oo pc, where ^200 is the distance 4 scaled 
to 200 pc). To produce the image, we subtracted the in- 
trinsic PN and particle background using a closed-filter 
observation in Small Window mode taken on 2002 Decem- 
ber 30. The raw data were binned in 3" x 3" pixels and 
smoothed with a Gaussian of FWHM = 15". Subtracting 
the background from a 1' x 1' box in the lower left corner 
of the PN image, we found the average image brightness of 
the tails (extraction box of 40" x 80" for each tail) about 

4 As mentioned by Kargaltsev et al. (2005), the Geminga's parallax reported by Caraveo et al. (1996) is incorrect, and the distance to 
Geminga is currently unknown; however, d ~ 200 pc seems to be a reasonable estimate. 
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agreement with the value reported by C03. 

The tails in the PN image look similar to those in the 
MOS images, reported by C03. Their shape fits the same 
analytical bow-shock model by Wilkin (1996), with an in- 
clination angle i w 90° between the line of sight and pul- 
sar velocity and a stand-off distance of about 20", hidden 
in the pulsar image broadened by the wide XMM-Newton 
PSF. We note, however, that although such one-zone mod- 
els can describe the FBS shape close to its apex, they are 
not expected to be applicable to PWN tails (see, e.g., B05). 

Thus, the PN data show the properties of the outer tails 
consistent with those derived by C03 from the MOS data. 
However, the patchy structure of the tails hints that there 
may be a substantial contribution from faint background 
objects (such as the point source 'X' resolved with Chan- 
dra; see Fig. 2). To separate such a contribution and prove 
that the tails are not an artifact, they should be observed 
with a better spatial resolution. 

3.2. Chandra ACIS Results 

The excellent resolution of Chandra provides a close- 
up view in the vicinity of Geminga (right panel of Fig. 
1) The image shows no emission 20" ahead of the pulsar, 
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Fig. 2. — ACIS-S3 image of a 1' X 5' region around Geminga, smoothed with a 7" FWHM Gaussian, with overlaid PN brightness contours 
(red curves). The polygon south of Geminga encloses the region of enhanced emission apparently connecting the arc (not seen at this resolu- 
tion) with the southern outer tail. The dashed rectangle shows the region used for the background evaluation. The 'X' marks a point source 
near the southern outer tail. 



predicted by C03. Instead, we see some diffuse emission 
at a distance of 5"-7", whose shape resembles an arc ex- 
tended perpendicular to the proper motion direction. The 
arc and the outer tails cannot be fitted with the same 
Wilkin's model. The arc-like polygon of a 116.5 arcsec 2 
area contains 32 counts, of which 10.1 counts are estimated 
to belong to the background (0.087 counts arcsec -2 , as 
measured in a 49" x 98" source-free rectangle north of the 
pulsar, shown in Fig. 2). This gives 21.9±5.7 background- 
subtracted counts. The arc's spectrum (T = 1.2 ± 0.4 at 
fixed TVjj = 1.1 x 10 20 cm" 2 , as obtained from a power- 
law fit using the C-statistic) is substantially harder than 
the pulsar's spectrum. A structure of such size and spec- 
trum cannot be ascribed to PSF tails of the pulsar's image, 
which does not show any significant pileup (the pulsar's 
count rate is only 0.07 counts/frame). The average image 
brightness of 1.0 x 10 -5 counts s -1 arcsec -2 corresponds 
to an intensity I arc ~ 0.9 x 10 -16 ergs cm -2 s -1 arcsec -2 , 
about 30 times brighter than the tails seen with EPIC. 
The X -ray luminosity of this structure is L arc ~ 5 x 10 28 

dloo er g s s_1 - 

The most striking feature in the ACIS image is a ~ 5"- 
wide axial tail, seen up to 25" (7.5 x 10 16 d2oo cm) from 
the pulsar in the direction opposite to the pulsar's proper 
motion. The tail is apparently detached from the pulsar 
by 5"-6"; its brightness is maximal at ~ 8", and it fades 
with increasing distance from the pulsar. The average im- 
age brightness in the 6" x 12" box (see Fig. 1, right panel), 
which contains 37.1 ± 6.6 background-subtracted counts, 
is 2.7 x 10 -5 counts s -1 arcsec -2 . Although the small 
number of counts precludes detailed spectral analysis, the 



tail's spectrum can be described by a power-law model 
with r = 1.0 ±0.2, somewhat harder than the spectrum of 
the outer tails. The intensity of the tail (in the 6" x 12" 
box), I ax iai = 3x 10 -16 ergs cm -2 s -1 arcsec -2 , is about 
2 orders of magnitude higher than that of the tails seen 
with EPIC. Its luminosity, L ax iai = (1.3 ± 0.2) x 10 29 d| 00 
ergs s -1 (measured from 49.1 ± 7.3 background-subtracted 
counts in a 6" x 20" box detached by 6" from the pulsar), is 
close to the total luminosity of the putative outer tails. An 
inspection of the EPIC PN image shows some patchy en- 
hancements along the tail direction, at a level of ^ 2 x 10 -7 
counts s -1 arcsec -2 , but they are indistinguishable from 
background fluctuations. 

Heavily smoothed ACIS images show an enhancement 
south of the pulsar which apparently connects the arc with 
the southern outer tail (an example is shown in Fig. 2). 
The region of enhanced emission within the 829 arcsec 2 
polygon includes 105 counts. Subtracting 71.9 counts of 
the scaled background, we obtain 33.1±10.8 excess counts, 
which corresponds to an intensity a factor of 6 higher than 
that in the EPIC tails. The slope of the excess spectrum, 
r = 1.0±0.4, is apparently similar to that of the axial tail. 
Thus, although the 3 a enhancement can hardly be con- 
sidered as a firm detection, its position and shape support 
the reality of both this structure and the southern outer 
tail. No enhancement is seen at the site of the northern 
outer tail; however, the northern tail looks less significant 
and more patchy in the EPIC images, and only a small 
part of it was imaged with ACIS. 
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4. DISCUSSION 

If the tails in the EPIC images are real, the Gcminga 
PWN is truly unique: a bow-shock-like structure with long 
outer tails and a short axial tail behind the pulsar have 
never been seen before in X-rays. Before discussing possi- 
ble interpretations of the observed structures, we note that 
the proper motion of Gcminga, 0'.'17 yr -1 , implies a pulsar 
speed v = 160 d km s _1 , where d = cfeoo/ sini. For a rea- 
sonable distance, it exceeds a typical sound speed in the 
interstellar medium (ISM), c s = 15 (/VO.6)" 1 / 2 ^ 172 km 
s" 1 , where ri and T = 10 4 T 4 K are the molecular weight 
and temperature. Assuming that the speed of a possible 
ISM flow at the location of Geminga is much lower than 
v, the ram pressure due to the pulsar motion in the ISM is 
Pram = 4.3 x 10~ 10 n d 2 ergs cm~ 3 , where n is the ISM den- 
sity in atomic mass units per cm 3 , This gives an estimate 
R h = 1.4 x 10 16 n _1 / 2 J _1 cm for the stand-off distance 
of the TS head, which translates into the projected an- 
gular distance IZh = 4'.'8 n _1 / 2 rf -2 . Thus, one can expect 
that Geminga is accompanied by a bow-shock PWN, with 
a characteristic size comparable to the sizes of the struc- 
tures observed with Chandra. We will discuss possible 
interpretations of the whole PWN, starting each from an 
assumption on the nature of the axial tail, the brightest 
feature of the PWN. 

4.1. The axial tail is a shocked pulsar wind? 

The axial tail could be interpreted as synchrotron emis- 
sion from the shocked PW collimated by the ram pres- 
sure. According to the simulations by B05, who assumed 
an isotropic PW, the TS has a bullet-like shape. For large 
Mach numbers, M — v/c s , and small values of the magne- 
tization parameter a of the pre-shock PW (see Kennel & 
Coroniti 1984), the bullet's cylindrical radius is rxs ~ -Rh, 
and the distance of its back surface from the pulsar is 
R h ~ 6R h . The shocked PW outside the TS is confined 
inside the CD surface which has a cylindrical shape be- 
hind the TS, with a radius tcd ~ 4i?h- The collimated 
PW flows with subrelativistic velocities: 0.1-0.3 c in the 
inner channel, r <; rxs, and up to 0.8-0.9 c in the outer 
channel, tts <5 t <^ tqd (see Figs. 1-3 and §3.3 in B05). 

First, one can speculate that the axial tail is the CD- 
confined cylindrical tube behind the TS, which implies a 
CD radius of <~ 3", IZh ~ 0'.'7sinz, and TZb ~ 5"sinz. In 
this interpretation, one should expect brightest emission 
from the shocked PW at <J 1" ahead of the pulsar, hidden 
within the pulsar image. The actually observed emission 
~ 5"-7" ahead of the pulsar (the arc) is not explained by 
this model. Being well outside the CD, the two outer tails 
cannot be associated with a shocked PW. One might spec- 
ulate that they are produced by the shocked ISM heated 
to X-ray temperatures, but the pulsar's speed is too low 
to support this speculation, and thermal models with rea- 
sonable temperatures do not fit the tails spectrum. Over- 
all, given the problems with explaining the observed PWN 
structure, this interpretation of the axial tail is hardly vi- 
able. 

Second, one could assume that the axial tail is asso- 
ciated with the shocked PW "sheath" immediately out- 
side the bullet-like TS, similar to the interpretations of 

5 Note, however, that the simulation of intensity by B05 do not si 



the tail behind PSR B1757-24 by Gvaramadze (2002) and 
the "tongue" of the X-ray Mouse PWN by Gaensler et al. 
(2004). The two outer tails might be associated with emis- 
sion from a shell immediately inside the CD surface, where 
the magnetic field is compressed (see B05) and hence the 
synchrotron emissivity could be enhanced 5 . In this inter- 
pretation, the outer tails should be parallel to the direction 
of pulsar motion, which is crudely in agreement with obser- 
vations, but their distance from the PWN axis, ~ 60"-70", 
implies a bullet diameter ~ 30"-35", much larger than the 
observed width of the axial tail. This interpretation also 
implies TZh ~ 16" sin i, which is larger than the distance 
to the observed arc unless the inclination angle is small, 
sinz <; 0.3. On the other hand, the width of the axial 
tail suggests TZh ~ 2'.' 5 sin i, smaller than the distance to 
the arc in the Chandra image. Moreover, the arc is, on 
average, a factor of 4 dimmer than the tail, and its extent 
perpendicular to the pulsar's proper motion is a factor of 4 
larger than the tail's width. Therefore, we are forced to as- 
sume that the TS head is unresolved from the pulsar, while 
the arc might be a CD head at a distance ~ 2.8i?h/ sini 
from the pulsar, larger than ~ 1.3i?h in the B05 simula- 
tions. The corresponding ratio rcr>/Rh ~ 28/ sini is much 
larger than ^4 in the B05 simulations. The simulations 
also do not explain the apparent detachment of the axial 
tail from the pulsar. These discrepancies could be caused 
by anisotropy of the pulsar outflow, neglected by B05. For 
instance, if the outflow is mostly equatorial (assuming the 
pulsar's spin axis aligned with its space velocity), the flat- 
tened TS head should be closer to the pulsar than in the 
isotropic case. 

We should also mention that this interpretation implies 
a rather large tail's length, l ta n ~ !)j ow r syn , where vn ow 
is the flow velocity and T S y n — 5.1 x lO^- 1 ^- 2 s is the 
synchrotron cooling time. The magnetic field just down- 
stream of the TS back boundary can be estimated, for 
a <C 1, as 

(cf. Kennel & Coroniti 1984). The synchrotron photons 
with maximum observed energies E w 8 keV are emitted 
by electrons with a Lorentz factor 

This gives 

/tail ~ 4X10 or c {—) 1 15*) ^7I Cm ' (3) 

which exceeds the observed length by two orders of mag- 
nitude even at lowest Ufl ow ~ 0.1c in the inner channel 
behind the TS, as found by B05. This discrepancy could 
be explained assuming that the surface brightness of the 
tail becomes too low at large distances from the pulsar 

bright shell-like structures. 
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(e.g., because of a decreasing magnetic field) to be seen in 
these images. 

Thus, the interpretation that the axial tail is the freshly 
shocked PW immediately outside the TS while the outer 
tails mark the CD surface is not quantitatively consistent 
with the available simulations. However, such an interpre- 
tation cannot be ruled out because those simulations do 
not take into account the intrinsic anisotropy of the PW. 

4.2. The axial tail is a pulsar jet? 

Another explanation of the axial tail, which we consider 
more plausible, is that it is a jet emanating from the pul- 
sar magnetosphere along the spin axis aligned with the 
pulsar's motion. The fact that only one jet is seen is not 
uncommon (PSR B1706— 44 is a vivid example; Ng & Ro- 
mani 2004), and it can be explained by Doppler boosting 
(the approaching jet is brighter than the receding one) 
and/or by intrinsic anisotropy of the polar outflows, or it 
may be caused by destruction of the forward jet by the 
ISM ram pressure. The Geminga's axial tail resembles the 
southeast jet ("inner counterjet" in Pavlov et al. 2003) of 
the Vela PWN, which is about twice brighter than the 
northwest jet in the direction of pulsar's motion. The 
Vela's southeast jet is also somewhat detached from the 
pulsar (perhaps because the polar outflow becomes visible 
only beyond a shock). Its projected length, 5 x 10 16 c?3oo 
cm, is close to that of the Geminga's axial tail, the spec- 
tra of both structures are very hard, and the ratios of their 
X-ray luminosities to the pulsar spindown powers are not 
very different: L x /E ~ 3.6 x 10 _6 <i| on for the Geminga's 
tail and ~ 0.6 x 10 — 6 c?§q for the Vela's southeast jet. 

If the jet is confined by its own magnetic field, a lower 
limit on the field can be estimated from the requirement 
that the electron Larmor radius is smaller than the jet 
radius, rj e t ~ 0.75 x 10 16 <i2oo cm, which gives 

b>90 (^v) 1/3 ^ /3aiG ' (4) 

where -Em is the maximum energy of the X-ray power-law 
spectrum (cf. Pavlov et al. 2003). It corresponds to the 
energy injection rate 

B 2 

W = g^(l + fc)vjlt 

where v- }e t — /3c is the bulk flow velocity in the jet, and 
k is the ratio of particle and magnetic energy densities. 
Even this lower limit on W is a substantial fraction of 
the spindown power for a mildly relativistic i>j e t expected 
(e.g., W ~ 0.05E for k = 1 and (3 = 0.5), which means 



that the magnetic field cannot strongly exceed the above 
lower limit. At such magnetic fields one would expect a 
jet length 

^jct ^jct^syn 

„ 1 r»i8 ft /100mG\ 3/2 /8keV\ l/2 

- 2 x 10 18 — — cm, 6 

0.5 V B J \ E M J ' K ' 

much larger than the observed ~ 7 x 10 16 <i cm. Therefore, 
we have to assume that the jet is destroyed or becomes 
uncollimated well before it loses its internal energy to ra- 
diation. 

If the axial tail is a pulsar jet, the outer tails could mark 
an equatorial outflow bent by the ram pressure. We are un- 
aware of theoretical PWN models which include both the 
PW anisotropy and the ram pressure effects. We, how- 
ever, expect that the equatorial PW component (which 
would produce a torus beyond a TS ring around a slowly- 
moving pulsar) would form a relatively thin shell between 
the TS and CD surface behind the pulsar, filled by a rela- 
tivistic plasma with a subrclativistic bulk flow velocity. If 
the outer tails turn out to be an artifact, we would sug- 
gest that most of the PW flows out of the magnetosphere 
along the spin axis while the equatorial PW component 
is unusually weak in Geminga (perhaps because of a small 
angle between the magnetic and spin axes). The arc ahead 
of the pulsar could be a head of the bent equatorial out- 
flow or remnants of a forward jet crushed be the ISM ram 
pressure. 

To summarize, the Chandra observation has conclu- 
sively shown the presence of PWN elements around the 
Geminga pulsar. With the sparse statistics of the Chandra 
data and the low spatial resolution of the XMM-Newton 
images, we cannot establish the nature of the extended 
emission unambiguously. However, it seems certain that 
the observed PWN structure implies that the Geminga's 
PW is intrinsically anisotropic. Much deeper Chandra ob- 
servations and modeling of magnetized anisotropic winds 
from fast-moving pulsars are needed to clarify the nature 
of this intriguing PWN 6 . 
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6 Aftcr our work was mostly completed, we became aware of the paper by De Luca ct al. (astro-ph/0511185) who reported on the same 
Chandra data. They found virtually the same properties of the axial tail but failed to notice the arc and the enhancement connecting the arc 
with the southern outer tail. 
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